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Abstract. In this paper we study the evolution of the two two-level àtoms interacting 
with a single-mode quantized radiation field, namely, two-atom multiphoton (, i.e. 
A;th-photon) Jaynes-Cummings model when the radiation field and àtoms are initially 
prepared in the superposition of displaced number states and excited atòmic states, 
respectively. For this system we investigate the atòmic inversion, Wigner function, 
phase distribution and entanglement. We show that there is a connection between all 
these quantities. Moreover, for symmetric (asymmetric) àtoms the system can generate 
asymmetric (symmetric) cat states based on the vàlues of the interaction parameters. 
This is shown in the behaviors of the Wigner function and phase distribution. The 
degree of entanglement for the field-atoms and the one-atom-remainder tangles depends 
on the energy foUow between the bipartite. Also it is sensitive to the interference in 
phase space and the value of the parameter k. 

PACS numbers: 42.50Dv,42.60.Gd 



1. Introduction 



Entanglement is the striking feature of quantum mechanics revealing the existence of 
nonlocal correlations among different parts of a quantum system. A pair of quantum 
systems is called entangled if the measurement on one of them cannot be performed 
independent of that of the other. Quantum entanglement is a resource for certain tasks 
that can be performed faster or in a more secure way than the classical correlation. For 
these reasons entanglement plays an essential role in quantum Information, e.g. quantum 
computing teleportation |0, cryptographic ||^, dense coding |^ and entanglement 
swapping . These new aspects have launched intensive experimental efforts to generate 
entangled states and theoretical efforts to understand their structures. For instance, the 
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entanglement between two qubits in an arbitrary pure state has been quantified by 
the concurrence 0, however, that of the mixed states has been given as the infimum 
of the average concurrence over all possible pure state ensemble decomposition. The 
closed form for the concurrence has been expressed by Peres- Horodecki measure 0]. 
Moreover, this measure has been extended to include a bipartite system AB, with 
arbitrary dimensions Da and Db in an overall pure state 0. The latter technique 
has been used to treat the entanglement for the tripartite quantum system in a Hilbert 
space with tensor product structure 2 2 oo P , e.g for the two-atom Tavis-Cummings 
model in an overall pure state. 

The interaction between the radiation field and the two-level atom, namely, Jaynes- 
Cummings model (JCM) is an important tòpic in quantum òptics since it is solvable 
in the framework of the rotating wave approximation (RWA) and it is experiment ally 
implemented UTO], p!ï[]. Also the JCM is a rich source for the nonclassical effects, e.g. the 



revival-collapse phenomenon (RCP) in the evolution of the atòmic inversion |jÏ2| and 
the generation of the cat states at one-half of the revival time ||r3|, |I4|. The importance 



of the JCM is increased as a result of the progress in the quantum Information |15 



The JCM has been generalized and extended in different directions |16|. One of these 



directions, which is of particular interest, is the two two-level àtoms interacting with a 
single quantized electromagnètic field (TJCM) |T^, [T^, ^ ^ . The atòmic inversion 
of the TJCM exhibits RCP having forms different from those of the JCM. This quantity 
has been investigated when the field is initially prepared in coherent state |T^, |T8|, p!9|, [20[], 
binomial state PT| and displaced squeezed state [E^]. Also the entanglement for the 



and binomial state [21 



TJCM when the field is initially prepared in coherent state 
have been discussed, too. 

The main object in quantum òptics is to increase the nonclassical effects obtained 
from the quantum system. This has been achieved by developing new states beside 



the Standard ones. One of these states is the displaced number state |]23[, which can be 
generated by passing coherent light through a non-linear medium [24|. In the framework 



of the superposition principle the superposition of displaced number states (SDN) has 
been developed as 



|a,m) = Ae[-D(a) + tD{—a)\\ni) 



E C(n,m)|n), 

n=0 



(1) 



where D{oí) is the displacement operator, e is a real parameter where its value will be 
specified in the text, is the normalization constant having the form: 



A: 



1 + + 2eexp(-2a2)L„(4a2 



(2) 
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and 

C{n,m) = X,{n\[D{a) + eD{-a)]\m). (3) 

The photon-number distribution P{n) related to (1) is 

P{n) = \C{n,m)\^ (4) 

For particular vàlues of tlie parameters the state (1) reduces to number states, coherent 
states, Schròdinger-cat states and displaced number state. The state (1) can be regarded 
as a single-mode vibration of electromagnètic field suddenly displaced by a coUection of 
two displacements vr out of phase with respect to each other. The generation of (1) has 
been estabhshed by the so-called "quantum state engineering" ||25|, Quite recently 



the decoherence for (1) in the framework of the Standard master equation, which is 
described by the phase insensitive attenuators or amplifiers, has been investigated in 
P7| . Also the superposition of the N displaced number states |^ and superposition 



of the squeezed displaced number states [|ï^, ^ have been developed, too. Finally, the 
density matrix p of the superimposed states has two parts, namely, statistical-mixture 
part Ps and interference part p/. The for mer provides Information on the original 
components of the states and the latter gives the interference in phase space. 

In the present paper we study the SDN against the two two-level àtoms interacting 
with the quantized single-mode electromagnètic field (TJCM). In the RWA the 
Hamiltonian [0, |T^, |T9|, which controls the system, takes the form 



H 



% = Hq + Hi, 
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(5) 



where Hq and Hj are the free and interaction parts of the Hamiltonian, àf and cr| are 
the Pauli spin operators of the jth atom; à {a^) is the annihilation (creation) operator 
denoting the cavity mode, u and Ua are the freqüències of the cavity mode and the 
atòmic systems (we consider that the two àtoms have the same frequency), \j is the 
atom-field coupling constant of the jth atom and k is the transition parameter. In 
the paper we mainly deal with the ratio g = A2/A1 and consider two cases, namely, 
symmetric and asymmetric according to g = 1 and g ^ 1, respectively. Furthermore, 
we assume that Ua = 2küü (, i.e. the exact resonance case) and the two àtoms and field 
are initially in the excited atòmic states |+, +) and the SDN, respectively. The atòmic 
ground state is denoted by |— ). Therefore, the dynamical state of the whole system can 



be evaluated as: 



|^(T)) = E C(n,m)[Xi(T,n,A;)|+,+,n)+X2(T,n,fc)|+,-,n + A;) 

n=0 



(6) 



+X3{T,n,k)\-,+,n + k) +X^{T,n, k)l 



n 



where T = Ait is the scaled time. The explícit forms for the dynamical coefíicients 
Xj{T, n, k) can be found in Jï^ |T^, ^ |2^, however, for analytical tasks we provide 
the forms of these coefíicients for g = 1: 



Xi{T,n, k) 



n\(n+k)\ 
[{n+fc)!]2+n!(n+2fc)! 



X2(r,n, fc) =X3(T,n, fc) 



(n+fc)! sm(TC„) 
n! Cn ' 



(7) 



X^{T, n, k) 



(n+fc)!y«!(«+2fc)! r . _ 

[(n+fc)!]2+n!(n+2fc)! [^^^1-^ SnJ -LJ , 



where 



[n + fc)! 



+ 2- 



2k)\ 



[n 



k)\ ■ 



(8) 



For this system we investigate the total atòmic inversion {az(T)) = |[(o"f(T)) + 
(cr|(T)), Wigner (W) function, phase distribution P(6) and entanglement between 
different components of the system. The object of this study is as foUows. There 
are some well-known phenomena occurred for the JCM so that we see the possible 
occurrence of such phenomena for the TJCM and what would be the influence of 
the interference in phase space on the behavior of the system. Also we quantify the 
entanglement between different components in the system. More illustratively, for JCM 
it has been shown that the Schròdinger-cat state can be generated at one-half of the 
revival time. Also this issue has been discussed only for the off-resonance TJCM when 
{k,g) = (1,1) with a » 1 in the framework of Q function [^]. Here we investigate 
this behavior for the resonance TJCM the most general case, i.e. for all vàlues of g, via 
W function, which gives better on quantum system than Q function. In this regard-for 
the system under consideration-we show that for the symmetric (asymmetric) case the 
asymmetric (symmetric) Schròdinger-cat states can be generated. This behavior likely 
occurs at the quarter of the revival time. Such generation of cat state is confirmed in 
the behavior of the P(0). Also we deduce the asymptotic form for the W function in 
the strong-intensity regime. Moreover, we shows that the degree of the entanglement is 
sensitive to the vàlues of both k, g, and e. Specifically, the degree of entanglement for 
the asymmetric case is greater than that of the symmetric one, and the interference in 
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phase space decreases the degree of entanglement. Finally, for k > 1 the evolution of 
the tangles can be connected by the behavior of the corresponding atòmic inversion. 

The paper is organized in the following order: In section 2 we investigate the 
evolution of the (^^(T)), W function and phase distribution. In section 3 we investigate 
the entanglement between different components of the system. In section 4 we 
summarize the main results. 



2. Atòmic inversion, Wigner function and phase distribution 

One of the bàsic quantities related to the JCM is the atòmic inversion, which represents 
the difference between the population of the excited and the ground atòmic states. The 
evolution of the atòmic inversion is representative by showing RCP. The occurrence of 
the RCP in this quantity has a quantum origin since it indicates the granular structure 
of the initial field distribution and also reveals the atom-field entanglement. A lot of 
efforts have been done for observing the RCP experimentally [IIÜ, O, RÜ, |3Ï], p2|. On 



the other hand, for the JCM it has been shown that there is a connection between the 
behaviors of the atòmic inversion and both the W function and phase distribution. In 
this section we discuss this situation for the TJCM and compare the obtained results 
with those of the JCM. Also we investigate the influence of the interference in phase 
space on the behavior of the system. 

We start with the total atòmic inversion, which can be evaluated for (6) as 
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Figure 1. (a) and (b) represent the (ctz(T)) against the scaled time T for {a,g) = 
(7,0.5), the curves A and B {{az{T)) + 2) in (a) and (b) are given for (e,77i) = 
(0, 0), (1, 0) and (0, 1), (0, 2), respectively. (c) P{n) against n for (e, a) = (0, 7) with the 
different vàlues of m as indicated and shifted from bottom by 0,0.1,0.2. The dashed 
curve in (c) represents the P{n) of the (e, a, m) — (1, 7, 0). 



oo 
n=0 



(9) 
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For the symmetric case and in the strong-intensity regime (SR), i.e. a » 1 and P{n) 
has a smooth envelope. The relation (y) reduces to 

i^ÀT)) = E |C(n,m)pcos (^^T^^^ . (10) 

It is evident that this formula is quite similar to that of the Standard JCM except the 
factor 3/2 has to be unity. Thus the behavior of the (^^(T)) of the TJCM is quite similar 
to that of the JCM, i.e. the atòmic inversion of the symmetric TJCM is insensitive of 
the existence of the other atom. In this case the revival time for (e, m) = (0,0), say, is 
Tr = 271^ n + 3/2. 

Now we draw the attention to the {cTzÍT)) of the asymmetric case (see Figs. 1 for 
given vàlues of the interaction parameters). Also we have plotted the corresponding 
initial P{n). From this figure one can see that when P{n) exhibits multipeak structure, 
e.g m 7^ 0, the revival patterns in the corresponding {az{T)) are split. This is related 
to that each peak of P{n) provides its own RCP in the evolution of the {az{T)), which 
interfere with those of the others to give such behavior. The comparison between the 
curves A and B in Fig. l(a) shows that the revival time of the initial cat states is 
approximately one-half of that of the initial coherent states and the locations of the 
symmetric and asymmetric revival patterns in the two curves are quite different. This 
can be easily realized by noting that each pj and ps has its own RCP. Furthermore, the 
revival time of the former is the two times smaller than that of the latter. Therefore, 
when the revival patterns contributed by pj and ps occur at the same interaction time 
asymmetric patterns are exhibited in {az{T)). This behavior is completely different 
from that of the symmetric case, which we have not presented. 

As is well known that for the JCM the evolution of the W function [ï^ as well as 



the Q function |3^, |3J, ^ provides Information on the RCP in the {aziT)). More 
illustratively, as the interaction is going on, the W function, which is initially presented 
by shifted Gaussian peak, splits into two peaks counter-rotate on a circle in the complex 
plane of the distribution with interference fringes in between. As the interaction proceeds 
the two peaks coUide at the opposite end of the circle having the initial form, they split 
again, and so forth. In the framework of the (cr^(T)): the col·lapse regions (revival 
patterns) in the Rabi oscillations occur in the course of the splitting (collision) of the 
distributions of the W function. In this regard, the Schròdinger-cat states are generated 



at the one-half of the revival time Also for the symmetric-off-resonance TJCM 

with the initial coherent light and using the Q function it has been numerically shown 
that asymmetric multipeak structure can be exhibited for a >> 1 Moreover, the 
widths and the heights of these peaks can be controlled by the vàlues of the detuning 
parameters and the type of the initial atòmic states. In this part we investigate the 
evolution of the W function for the TJCM when a is relatively small. This is inspired 
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by the facts that the W function gives information on the quantum system better than 
Q function and the cat states generated for small a are more nonclassical than those 
generated for strong initial intensity [^]. Additionally, we will deduce the asymptotic 
from for the state generated by the system when a » 1. The W function for (6) can 



be easily evaluated as |I4 



n,n'=0 

'{2\z\') + (-l)nX2(n,T)X2(n',r) +X3(n,r)X3K,T)]y^L— '(2|^n (11) 



+X,in,T)X,{n',T)J^^L:7,Í,i2\z\')}, 



where z = x + iy and L'^i·) is the associated Laguerre polynomial. Generally, we have 
noted for small vàlues of a that the entanglement between the two àtoms through the 
bosonic system leads to that the initial form of the W function cannot appear again 
after switching on the interaction. For instance, for the symmetric case the single-peak 
structure coming from the leading terms (cf. (7)) is dominant with distortion around 
circle in the xy-plane (see Figs. 2). In Figs. 2 the vàlues of the interaction time have 
been chosen from the evolution of the corresponding (o"^(T)). Comparison between Figs. 
2(a), (b) and (c) shows that the nonclassicality (, i.e. the negativity) is more pronounced 
when T = Tr/A. Additionally, the negative vàlues in W function at T = T,. is much 
greater than those at T = Tr/2, however, subsidiary peaks start to appear for the latter 
case. This will be remarkable in the behavior of the P{Q), as we will see below. 

Now we derive the asymptotic form for the W function for the case a >> 1 and 
{g, k, e, m) = (1, 1, 0, 0) using SR. In this case we can replace C(n, 0) by C{n+ 1, 0) and 
C(n+2, 0) for the second and third terms in the curely brackets in (11), respectively, and 
then apply the SR for the arguments of the trigonometric functions. After some tricks. 
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using the gener ation function of the Laguerre polynomial and the Taylor expansion, the 
expression (11) can be evaluated as 




Figure 3. The W function against x and y for the asymptotic form (12) {a,T) = 
{7,Tr/A) (12) (a); the exact form (11) for {g,a,e,T) = (0.5,3,0,19.47003) (b) and 
(0.5,3,0,10.18) (c). 



W{x,y,T) ~ ^|2exp[— í/^ — {x — V^a^] + exp[— (x — v^acosr;)^ — {v — V^asinr/)^] 

(12) 

+ exp[-(a; - \/2acos?7)^ - {v + v^asinr;)^] - 2Iint{T) sin r^j, 



where rj = the interference part hntiT) takes the form 

Iint{T) = exp[— (x — \/2a cosr])"^ — y'^] cos[2?7' + 2\/2xa sinr] — Q;^sin(2í;)] sinr; 

+ exp[-(x - V^acos^ f )^]{exp[-(y + ^smrjf] sm//+ (13) 

exp[-(y - -^sinriy] sin//_| 
and 

T 

Ijl± — rj' — sin rj + V2a(x sin 77 ± y cos 77) ± V2ay, ri' = TV^ + (14) 

The formula (12) can give Information on the W function for large vàlues of a for 
which it is diíïicult to deal with the exact form (11) because of the associated Laguerre 
polynomial. Prom (12) it is obvious that the first term is time independent since it 
is connecting with the leading terms in (7). Also W function includes three individual 
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peaks and three interference fringes each of them locating between two of these individual 
peaks. When a is very large T,^ ~ 2-n\fE and we can extract various facts from 
(12). For instance, the interference part Iint{T) has no (maximum) contribution at 
T = Tr,Tr/2 (Tr/4). To be more specific, when T = Tj. the W function reduces to 
that of the initial one, which is Gaussian bell centered at (v^a, 0), whereas at T = Tj./2 
the initial bell sphts into two typical peaks localized at (±\/2qí, 0) and eventually at 
T = Tr/A three-peak structure centered at (A/2a, 0), (0, ±\/2a) with interference fringes 
in between are exhibited, i.e. three-component cat state is generated (see Fig. 3(a)). 
Moreover, we have checked the Q function for these vàlues of the interaction times 
with a >> 1 and obtained number of peaks as those exhibited for the W function. 
This behavior of the W function is different from that of the JCM even though the 
evolution of the atòmic inversion for the two systems is quite similar. This can be easily 
understood by comparing the generalized Rabi oscillation for the JCM and the TJCM. 
Actually, for the symmetric case the Rabi oscillation of the TJCM is approximately two 
times greater than that of the JCM (cf. (|])). This fact explains why the cat states in 
the TJCM are generated at the quarter of the revival time, however, in the JCM at the 
one-half of the revival time. On the other hand, we have plotted the W function (11) for 
the asymmetric case in Figs. 3(b) and (c) for the revival and col·lapse times, respectively, 
when a = 3. From these figures one can observe that the W function exhibits symmetric 
shapes in phase space, which are in contrast with the the symmetric case (see Figs. 2). 
The nonclassical effects are more pronounced in the course of the revival time, which 
also are greater than those occurred for the symmetric case (compare Figs. 3(b) and 
(c) as well as Figs. 3 (b)-(c) and Figs. 2). Furthermore, the shapes presented in Figs. 
3(b)-(c) indicate that the generated cat states are of the microscopic type. The W 
function of the microscopic cat state exhibits a complicated shape resulting from that 
the contribution of the different components of the state are located close to the phase- 
space origin. It is worth mentioning that the microscopic cat state (, i.e. a is small) 
provi des more nonclassical effects than the macroscopic one (compare Figs. 3(a) 
and (b)). 

Now we turn the attention to the phase distribution P(9). To investigate the P(9) 
we use the Pegg-Barnett formulism |3^, which for the density matrix elements p„^„/(T) 
gives 
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Figure 4. The phasc distribution P(0) against 8. Figure (a) is given for 
(a,e, fc) = (3,0,1) with (T,.g) = (T^,!) (sòlid curvc), (r^/2, 1) (short-dashcd curve), 
(T^/4,1) (long-dashed curve), (19.47003,0.5) (star-centered curve) and (10.18,0.5) 
(bell-centered curve). For (b) {a,e,m,g,k) = (7,1,0,0.5,1) with T = 40.99995 
(second-revival time, sòlid curve), 9.099998 (coUapse-time, short-dashed curve) and 
21.00004 (first-revival-timc, long-dashed curve). In (b) the star-centered curve is given 
for {a,e,m,g,k,T) = (7,0.1,0.5,1,41.04966). For (c) {a,e,m,g,k) = (7,1,0,1,2) for 
T = 7r/4 (sòlid curve), n/2 (short-dashed curve) and tt (long-dashed curve). In (c) 
star-centered curve is given for (a, e, m, g, k, T) = (7, 0, 0, 1, 3, 7r/4). 



n,n'=0 



= è|l + 2 E C{n,rn)C{rí,m)[X^{n,T)X^{rí,T) + X2{n,T)X2{rí,T) (15) 

^ n>n' 



+X3(n, T)X^{n', T) + ^4(71, T)Xi{n' , T)\ cos[(n - n')e]}, 



where we have considered that the phase reference ©0 = 0. It is evident that 
P(©) = P(— ©). As a result of this symmetry we have plotted P(©) against © for 
< © < TT in Figs. 4(a)-(c) for given vàlues of the interaction parameters. Prom Fig. 
4(a), i.c. for the symmetric case, one can see that in the course of the revival time the 
P(0) exhibits single-peak structure at © = 0. This behavior is close to that of the 
initial coherent hght. Nevertheless, at one-half (quarter) of the revival time two small 
lateral peaks (wings) occur. This indicates generation of the asymmetric cat states in 
the system. This behavior is not clear in the corresponding W function (see Figs. 2). It 
seems that the non-occurrence of the lateral peaks in the W function is related to that 
a is considerably small. In this case the contribution of the diíferent components in the 
W function will be close to the phase-space origin and destructively interfere with each 
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others. Furthermore, from the star-centered and bell-centered curves in Fig. 4(a) (, 
i.e. asymmetric case) P(G) exhibits two-peak (four-peak) structure through the revival 
(coUapse) time. In other words, two-component and four-component cat states can be 
generated in the asymmetric TJCM based on the vàlues of the interaction parameters. 
This has been explained for the W function as a generation of microscopic cat states. 
Information on the interference in phase space is shown in Fig. 4(b) for the asymmetric 
case. For the asymmetric case there are two forms of the revival patterns: one is smooth 
and the other is asymmetric. Thus in Fig. 4(b) we have plotted the P(G) for the first 
and second revival times. As is well known that for Oq = the -P(O) of the cat state 
possesses single-peak at O = and two wings as ^ ±7r This form changes by 
switching on the interaction where in the course of the collapse time asymmetric eight 
peaks are created, which are transformed to symmetric four peaks at the first revival 
time and eventually are collapsed to two-peak structure at the second revival time. It 
is obvious that the shape of the peaks in P(0) includes Information on the shape of 
the revival patterns, i.e. when the P(0) exhibits smooth peaks the revival patterns are 
smooth and vice versa. Furthermore, the deformation in the P(n) leads to deformation 
in the peaks of P(0) (see the star-centered curve in Fig. 4(b)). Also we have noted 
that when the vàlues of the ratio g changes the number of peaks in P(0) changes, too. 
Furthermore, we have noted that the behavior of the case = 2 is similar to that of 
k = 1, however, for the former the peaks are always smooth and narrower than those for 
the latter (see Fig. 4(c)). Also when k > 2 the P(0) exhibits multipeak structure (see 
the star-centered curve in Fig. 4(c)). This situation is similar to that of the two- mode 



JCM g|. 



3. Entanglement 

Entanglement is at the heart of quantum Information theory. Various efforts are 
done to characterize qualitatively and quantitatively the entanglement properties of 
the quantum systems. This has been motivated by the progress in the experimental 
techniques aiming to create entangled states, which provide phenomena very different 
from the classical physics In this section we investigate the the entanglement for 



the system under consideration using the tangle / defined in ||^ . Based on the symmetry 
exchange for the system we investigate two forms of tangles, which are field-atoms and 
one-atom-remainder. In other words, assume that /, Ai and A2 denoting the field, first 
atom and second atom, respectively. Therefore, the field-atoms tangle (/ — A1A2) and 
one-atom-remainder tangle {Ai — /A2), say, are define as 0: 

If-A.A,{T) = 2[1 - Trp^(T)] = 2[1 - Trpi^^^(T)], 

(16) 

lM-fMÍT) = 2[1 - Trpi^(T)] = 2[1 - Trp2^^(T)], 



12 



where, e.g., Pf{T) is the reduced density matrix of the field, which can be obtained 
by tracing the total density matrix of the system over a complete set of the àtoms 
A1A2. The forms (16) quantify the degree of entanglement to which the ensembles 
behave as a coUective entity. It is worth mentioning that the tangle here includes the 
notion of the purity. Generally, when If^AiA2{T) = 0, say, the parties / and A1A2 are 
completely disentangled, but of course could be in states difïerent from those of the 
initial ones. Nevertheless, when /j_yi^^2(T) = 2 the parties are maximally entangled. 
In the foUowing we use the terminologies: If_AiA2ÍT) < 1 for weak entangled parties 
and If-AiA2{T) > 1 for strong entangled parties. 
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Figure 5. The tangle I f-AiA-iiT) against the scalcd time T when a = 7 for (a) 
(e, g, m, k) = (0. 1, 0, 1) (the short-dashed curve), (0, 0.5, 0, 1) (the long-dashed curve), 
(1,1,0,1) (the solid-curve-1) and (0,0.5,2,1) (the sohd-curve 2); (b) {€,g,m,k) = 
(0,1,0,2) (short-dashed curve) and (0,0.5,0,2) (sohd curve); and (c) {e,g,m) = 
(0, 1, 0). The solid-curve-2 in (a) is shifted by 1. 

3.1. Field-atoms tangle 

In this part we study the dynamical evolution for the tangle If^AiA2{T) for the system 
under consideration, which can be easily evaluated as 

If_A,A2{T) = 2-2 E \C{n, m)C{n', m)Xi{n, T)Xi{n', T) + C{n - k, m)C{n' - k, m) 

n,n'=0 ^ 

X [X2Ín - k, T)X2{n' - k,T) + Xs{n - k, T)X^{n' - k, T)] (17) 
+C(n - 2k, m)C(n' - 2k, m)X4(n - 2k, T)Xi{rí - 2k, T)}^ . 

We have plotted (17) in Figs. 5(a)-(c) for different vàlues of the interaction parameters. 
For all curves one can observe that at T = 0, where the bipartite is disentangled, 



13 



If^AiA2{T) = 0. We start with the symmetric case in Fig. 5(a). From the short- 
dashed curve it is obvious that when í > the entanglement is established between 
the two parties and If-AiA2(T) goes rapidly to the strong-entangled region, stays for 
while, dips to weak entangled region and eventually becomes steady for large interaction 
time. This behavior is completely different from that of the purity of the JCM, which 
exhibits oscillatory behavior in a good correspondence with the revival patterns in the 



atòmic inversion, and reduces to a pure state at the middle of the collapse time fTi 
From the sohd-curve-1 one can observe that the interference in phase space decreases 
the degree of entanglement at particular vàlues of the interaction time and increases the 
oscillatory behavior in the tangle (compare the short-dashed curve and sòlid- curve-1). 
For the asymmetric case, i.e. the long-dashed curve and solid-curve-2, one can observe 
generally that the degree of entanglement becomes greater than that of the symmetric 
case (compare short-dashed and long-dashed curves). Also the influence of including 
Fock state in the optical cavity on the behavior of the If-AiA2{T) is shown by the 
solid-curve-2. From this curve it is obvious that the oscillatory behavior in /^.^^^^(T) 
is increased compared to that of the initial coherent light, however, this behavior is 
not consistent with the RCP in the corresponding atòmic inversion (see Fig. l(b)). 
Now we draw the attention to Fig. 5(b) in which we have plotted If-AiA2(T) for the 
two-photon transition case, i.e. k = 2. This figure indicates that the bipartite can be 
periodically disentangled for the symmetric and asymmetric cases. For the symmetric 
case the If^AiA2(T) stays in the strong-entangled region most of the interaction time 
and goes to zero periodically with period vr. This can be easily understood, e.g. for 
{e,g,m) = (0, 1,0), in the framework of the SR as follows. In this regime such terms 
1/n, 1/n^, .. tend to zeros leading to that Cn — 2(n + 5/2) (cf. (||)). Thus the state (6) 
can be expressed in the following asymptotic form: 



l^(r)) - i E C(n,0){[cos(2nT + 5T) + 1] |+, +, n) - 2i sin(2nT + 5T) 

(18) 

x[|+,-,n + 4) + |-,+,n + 2)] + [cos(2nT + 5T) - 1]|-, -, n + 4)}. 



When T = stt, s is odd number, the state (18) reduces to 



|^(T)) ^ C'í^' 0)1^ + 4) ® exp(m)|-, -). (19) 

n=0 



It is evident that in ( |T9| ) the state of the radiation field is quite similar to that of the 
initial one but with four-photon shift. This behavior is similar to that of the purity 
of the JCM, but the forms of the phase factors and the number of shifted photons in 
the two quantities are different (compare (|T^ with (22)-(23) in [0]). Moreover, when 



T = s'tt where s' is even integer, (18) reduces to its initial form. Also at T = vrs"/2, s" 
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Figure 6. The /ai-/A2(^) against the scaled time T when a = 7 for (a) (e, g, m, k) = 

(0, 0.5, 0, 1) (curve A), (0, 1, 0, 1) (curve B), (0, 0.5, 2, 1) (curve C) and (1, 1, 0, 1) (curve 
D), and (b) {e,g,m,k) = (0,1,0,2) (sòlid curve) and (0,0.5,0,2) (dashed-curve). In 
(a) the curves A-D are shifted from bottom by 0, 0.3, 0.7, 1.4, respectively. 

is odd integer, the state (18) takes the form 

\^{T)) 1 E C{n,0){[\+,+,n) - |-,-,n + 4) 

n=0 , , 

(20) 

-2í(-l)«[|+,-,n + 2) + |-,+,n + 2)]}. 

For (20) one can easily prové that If^AiAziT) — 1- AU these analytical facts are 
remarkablc in Fig. 5(b). For the asymmetric casc, i.c. the sohd curve in Fig. 5(b), 
we obtain behavior as that of the symmetric casc but the degrec of entanglement is 
grcatcr and the disentanglement pcriod is 27i. Nevertheless, for this casc it is difficuh 
to obtain cxphcit forms for the statc vector at If-AiA2{T) = 0. It is obvious that 
the disentanglement period depends on the value of the ratio g. Also the behavior 
of the If-AiAiiT) is consistent with the occurrence of the RCP in the evolution of 
the corresponding ((T^(T)), we have checked this fact. From Fig. 5(c) one can see 
that for A; = 3 the bipartite are maximally entangled, however, for A; = 4 the degree 
of entanglement is decreased and the bipartite becomes periodically disentangled with 
period tt. This behavior is similar that of the purity of the JCM. Finally, in Fig. 5(c) 
we have plotted If-AiA2 (T) of the asymmetric case only where that of the symmetric 
case provides quite similar behaviors. 
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3.2. One-atom-remainder tangle 

In this part we study the dynamical evolution for the tangle /^^_/^2(T), which can be 
expressed as: 

iM^fA-AT) = 2-2 E {|C(n,m)CKm)n|Xi(n,T)|2 + |X3(n,T)p][|Xi(n',T)|2 + |X3K, 

n,n'=0 ^ 

+ |C(n,m)C(n',m)|2[|X2(n,T)|2 + |X4(n,T)|2][|X2K,T)|2 + |X4KT)|2] 

+2C{n + k, m)C{n, m)C{n' + k, m)C{n', m) [Xi(n + k, T)X;{n, T) + X^n + fc, T)XAn, T)] 

x[Xi{rí + k,T)X*{rí ,T) + X^rí + A;,T)X4(n',T)]}. 

We have plotted lA^-fA2 (T) in Figs. 6 for the given vàlues of the interaction parameters. 
The obvious remark from these figures is that < lAi~fA2{T) < 1- This is completely 
different from the behavior of the If^AiA2{T) (compare Figs. 5 and 6 as well as 
expressions (17) and (21)). This indicates that the degree of entanglement for the one- 
atom-remainder tangle is less than that of the field-atoms tangle. Also this shows that 
the quantum entanglement cannot be equally distributed among many different objects. 
This can be explained as foUows. We have an isolated system, i.e. the interaction with 
the environment is neglected, and hence the energy in the field and the two-atom system 
is periodically exchanged. In other words, when the k photons are annihilated from the 
radiation field they are created, i.e. equally distributed, in the two-atom party and 
vice versa. This means that the energy included in the fA2 party is more than that in 
the Al party. Actually, the entanglement is a direct consequence of the energy foUow 
between the two parties. The rate of foUow of energy in the bipartite of the If-AiA2(T) 
is greater than that in the lAi-fA2ÍT) and this could be the origin in the difference 
between the evolution of the two tangles. Now from the curves B and D in Fig. 6(a), 
i.e. symmetric case, one can observe that in contrast to If-AiA2(T) the lAi-fA2(T) 
exhibits an oscillatory behavior in a good correspondence with the revival patterns in 
{az{T)). Also the interference in phase space decreases the degree of entanglement in the 
lAi~fA2{T) (compare curves B and D). It is worth mentioning that the curve D is quite 
similar to that of the purity of the JCM (see Fig. 2 in [0). For the asymmetric case, and 
from the curve A one can observe that lAi~fA2{T) exhibits compound behavior, i.e. it 
includes periòdic and oscillatory behaviors. Also the bipartite can be approximately 
disentangled in the course of the first revival pattern. As in If_AiA2ÍT) the Fock 
state increases the range of the oscillatory behavior in the lAi-fA2{T) (see the curve 
C). From Fig. 6(b), which is given for the case k = 2, one can see that for the 
symmetric case the lAi-fA2ÍT) exhibits disentanglement periodically with period vr (cf. 
(|I^)) as in If_AiA2{T) and also includes oscillatory behavior. The lAi-fA2ÍT) of the 
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asymmetric case (, i.e. dashed curve) provides two types of disentanglement long lived 
(at T — 7r/2,37r/2) and instantaneous (at T — 27r) and they occur periodically with 
period 2tt. For the cases /c = 3, 4 we have noted that the lAi-fA2{T) — 1 with periòdic 
disentanglement only for the case k — A. 

4. Conclusion 

In this paper we have trcatcd the system of two two-level àtoms intcracting with 
multiphoton single-mode field. The two àtoms and the field are initially prepared 
in the excited atòmic states and in the superposition of displaced number states, 
respectively. We have investigated the behavior of the atòmic inversion, W function, 
phase distribution and entanglement. We have shown that the TJCM can generate 
asymmetric (symmctric) cat states at quarter of the revival time for the symmetric 
(asymmetric) casc. Also this has been confirmed in the behavior of the phase 
distribution. Moreover, we have deduced the asymptotic form for the W function when 
Q! >> 1. Also we have noted that when the vàlues of g changes the number of the 
components of the generated cat state in the system changes, too. The P(6) exhibits 
multipeak structure íor k > 2. For the entanglement we have investigated two types, 
namely, field-atoms and one-atom-remainder tangles. We have shown that the degree of 
entanglement in If_AiA2 is much greater than that in lAi-fA2- We have explained this 
in the framework of the rate of energy follow between differcnt parties of the system. We 
have obtained the foUowing facts related to entanglement. The degree of entanglement 
for the asymmetric case is greater than that of the symmetric one. The interference 
in phase space decreases the degree of entanglement. There is a similarity between the 
dynamical behavior of the If_AiA2{T) and the purity of the JCM only for /c > 1. 

Acknowledgments 

I would like to thank the Abdus Salam International Centre for Theoretical Physiscs, 
Strada Costicrs, 11 34014 Trieste Italy for the hospitality and financial support under 
the system of associateship, where a part of this work is done. 

References 

[1] Benenti G, Casati G and Strini G 2005 "Principle of Quantum Computation and Information" 

(World Scientific, Singapore). 
[2] Bennct C H, Brassard G, Crepeau C, Jozsa R, Peresand A and Wootters W K 1993 Phys. Rev. 

Lett. 70 1895. 

[3] Ekert A 1991 Phys. Rev. Lett. 67 661; Cirac J I and Gisin N 1997 Phys. Lett. A 229 1; Fuchs C 
A, Gisin N, Griffiths R B, Niu C-S and Peres A 1997 Phys. Rev. A 56 1163. 



17 



[4] Ye L and Guo G-C 2005 Phy. Rev. A 71 034304; Mozes S, Oppenheim J and Reznik B 2005 
Phys. Rev. A 71 012311. 

[5] Glòckl O, Lorenz S, Marquardt C, Heersink J, Brownnutt M, Silberhorn C, Pan Q, Loock P V, 
Korolkova N and Leuchs G 2003 Phys. Rev. A 68 012319; Yang M, Song W and Cao Z-L 2005 
Phys. Rev. A 71 034312; Li H-R, Li F-L, Yang Y and Zhang Q 2005 Phys. Rev. A 71 022314. 

[6] Wootters W K 1998 Phys. Rev. Lett. 80 2245; Hill S and Wootters W K 1997 Phys. Rev. Lett. 
78 5022. 

[7] Peres A 1996 Phys. Rev. Lett. 77 1413; Horodecki P 1997 Phys. Lett. A 232 333. 
[8] Rungta P, Buzek V, Gaves G M, Hillery H and Milburn G J 2001 Phys. Rev. A 64 042315. 
[9] Tessier T E, Deutsch I H and Delgado A 2003 Phys. Rev. A 68 062316. 
[10] Kaluzny Y, Goy P, Gross M, Raimond J M and Harochc S 1983 Phys. Rev. Lett. 51 1175; Rempe 

G, Walthcr H and Klein N 1987 Phys. Rev. Lett. 57 353. 
[11] Boca A, Miller R, Birnbaum K M, Boozer A D, McKeever J and Kimble H J 2004 Phys. Rev. 
Lett. 93 233603. 

[12] Eberly J H, Narozhny N B and Sanchez-Mondragon J J 1980 Phys. Rev. Lett. 44 1323; Narozhny 

N B, Sanchez-Mondragon J J and Eberly J H 1981 Phys. Rev. A 23 236; Yoo H I, Sanchez- 
Mondragon J J and Eberly J H 1981 J. Phys. A 14 1383; Yoo H 1 and Eberly J H 1981 Phys. 
Rep. 118 239. 

[13] Gea-Banacloche J 1991 Phys. Rev. A 44 5913; Buzek V, Moya-Cessa H and Knight P L 1992 
Phys. Rev. A 45 8190; Buzek V and Hladky 1993 J. Mod. Opt. 40 1309; Zaheer K and Wahiddin 
M R B 1994 J. Mod. Opt. 41 151. 

[14] Faisal A A El-Orany and Obada A-S 2003 J. Opt. B: Quant. Simeclass. Opt. 5 60. 

[15] Bose S, Fuentes-Guridi I, Knight P L and Vedral V 2001 Phys. Rev. Lett. 87 050401; Ibid 279901. 

[16] Barnett S M and Knight P L 1986 Phys. Rev. A 22 2444; Sharma M P, Cardimona D A and 
Gavriclidcs A 1989 Opt. Commun. 72 291; Cardimona D A 1990 Phys. Rev. A 41 5016; Abdalla 
M S, Ahmed M M A and Obada A-S F 1990 Physica A 162 215, and references therein. 

[17] Iqbal M S, Mahmood S, Razmi M S K and Zubairy M S 1988 J. Opt. Soc. Am. B 5 1312. 

[18] Sharma M P, Cardimona D A and Gavrielides A 1989 J. Opt. Soc. Am. B 6 1942. 

[19] Jcx I, Matsuoko M and Koashi M 1993 Quant. Opt. 5 275 

[20] .Jcx I 1990 Quant. Opt. 2 443 

[21] Abdel-Aty M 2004 J. Opt. B: Quant. Semiclass. Opt. 6 201. 
[22] Xu L, Zhang Z-M and Ghai J-L 1991 J. Opt. Soc. Am. B 8 1157. 
[23] Olivcira F A M, Kim M S, Knight P L and Buzek V 1990 Phys. Rev. A 41 2645. 
[24] Olivcira G C D, Almeida A R D, Queirós I P D, Moraes A M and Dantas C M A 2005 Phys. A 
351 251. 

[25] Moya-Cessa H 1995 J. Mod. Opt. 42 1741; Dantas C M A, Queiroz J R and Baseia B 1998 J. 
Mod. Opt. 45 1085. 

[26] Vogel K, Akuhn V M and Schleich W P 1993 Phys. Rev. Lett.71 1816; Parkins A S, Marte P, 

ZoUer P, Carnal O and Kimble H J 1995 Phys. Rev. A 51 1578. 
[27] Dodonov V V and L A de Souza 2005 J. Opt. B: Quant. Semiclass. Opt. 7 S490. 
[28] Obada A-S F and Abd Al-Kader G M 1999 J. Mod. Opt. 46 263; Faisal A A El-Orany 1999 Czeh. 

J. Phys. B 49 1145; Faisal A A El-Orany Pcfina .J and Abdalla M S 1999 J. Mod. Opt. 46 

1621; Kamli A and Bougouffa S 2005 Act. Phys. Hung. B 23 1. 
[29] Kim M S and Agarwal G S 1999 J. Mod. Opt. 46 2111. 



18 

[30] Meunier T, Gleyzes S, Maioli P, Auffeves A, Nogués G, Brune M, Raimond J M and Haroche S 

2005 Phys. Rev. Lett. 94 010401. 
[31] Faisal A A El-Orany 2004 ./. Phys. A: Math. Gen. 37 6157; Faisal A A El-Orany 2005 J. Phys. 

A: Math. Gen. 38 5557; Faisal A A El-Orany 2005 J. Opt. B: Quant. Semiclass. Opt. 6 432; 

Faisal A A El-Orany 2005 J. Mod. Opt. (submitted); Faisal A A El-Orany 2005 J. Phys. B: 

At, Mol, Opt. (submitted); Faisal A A El-Orany 2005 Opt. Commun. (submitted). 
[32] Faisal A A El-Orany 2004 J. Phys. A: Math. Gen. 37 9037. 
[33] Eiselt J and Risken H 1989 Opt. Gommun. 72 351. 
[34] Eiselt J and Risken H 1991 Phys. Rev. A 43 346. 

[35] Eiselt J and Risken H 1991 Phys. Rev. A 44 4623; Werner M J and Risken H 1991 Quant. Opt. 
3 185. 

[36] Millcr C A, Hilscnbcck J and Risken H 1992 Phys. Rev. A 46 4323. 
[37] Faisal A A El-Orany 2002 Phys. Rev. A 65 043814. 

[38] Faisal A A El-Orany, Pefina J and Abdalla M S 2000 J. Opt. B: Quant. Semiclass. Opt. 2 545. 
[39] Pegg D T and Barnett S M 1988 Europhys. Lett. 6 483; Ibid 1989 Phys. Rev. A 39 1665; Barnett 

S M and Pegg D T 1989 J. Mod. Opt. 36 7. 
[40] Faisal A A El-Orany, Mahran M H, Wahiddin M R B and Hashim A M 2004 Opt. Gommun. 240 

169. 

[41] Bouwmeester D, Ekert A and Zeillinger A 2000 "The physics of quantum Information" (Springer, 
Berlin) . 

[42] Rekdal P K, Skagerstam B-S K and Knight P L 2004 J. Mod. Opt. 51 75. 



